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Abstract 



(-H ' The heat operator with a pure soliton potential is considered and its Green's func- 

tion, depending on a complex spectral parameter k, is derived. Its boundedness prop- 
erties in all variables and its singularities in the spectral parameter k are studied. A 
generalization of the Green's function, the extended resolvent, is also given. 

^2 ' 1 Introduction 



The Kadomtsev-Petviashvili equation in its version called KPII 

(u t - 6uu Xl + u XlXlXl ) Xl = -3u X2X2 , (1.1) 



where u = u(x,t), x = (xijXa) and subscripts x\, xi and t denote partial derivatives, is 
a (2+l)-dimcnsional generalization of the celebrated Kortcwcg do Vrics (KdV) equation. 
The KPII equations, originally derived as a model for small-amplitude, long- wavelength, 
weakly two-dimensional waves in a weakly dispersive medium pQ, have been known to be 
integrable since the beginning of the 1970s [2] [3], and can be considered as a prototypical 
(2+l)-dimensional integrable equation. 

The KPII equation is integrable via its association with the operator 

C(x,d x ) = -d X2 +dl 1 -u(x), (1.2) 

which define the well known equation of heat conduction, or heat equation for short. The 
spectral theory of the operator (|1.2|) was developed in [l]-[2] in the case of a real potential 
u(x) rapidly decaying at spatial infinity, which, however, is not the most interesting case, 
since the KPII equation was just proposed in [1] in order to deal with two dimensional weak 
transverse perturbation of the one soliton solution of the KdV. 

A spectral theory of the KPII equation that also includes solitons has to be build using the 
resolvent approach, as already successfully done for the KPI equation [8; . In this framework 
it was possible to develop the inverse scattering transform for a solution describing one 
soliton on a generic background [5J , and to study the existence of the (extended) resolvent 
for (some) multisoliton solutions [ID] . However, the general case of A-solitons is still open. 



1 



Following [5], the first step in building the inverse scattering theory for this case lies in 
building a Green's function Q(x,x','k) of the heat operator (11.21) corresponding to the pure 
soliton potential u(x), that is, satisfying 



(~d X2 + dl - u{x)) G(x, x', k) = 5(x - x'), (1.3) 

and such that 

G(x,x',k) = e^^i^'J+k 2 ^-^) g^^k) (1.4) 

is bounded with respect to the variables x, x' G R 2 and k 6 C and has finite limits at infinity. 
The construction of such function is the subject of this article. 

2 Heat operator with pure soliton potential and its Jost 
solutions 

Soliton potentials (see [l0]-[16] for details) are labeled by the two numbers (topological 
charges) N a and Nb, which obey condition 

N a ,N b >l. (2.1) 

Let 

N = N a + N b , (2.2) 
so that AT > 2. We introduce the TV" real parameters 

Kl < K2 < ■ ■ ■ < km, (2.3) 

and the functions 

K n (x) = n n x 1 + n 2 n x 2 , n = l,...,M. (2.4) 

Let 

e K( -^ = di&g{e K ^}% =1 (2.5) 

be a diagonal TV" xW" matrix, let V be a TV x iVf, constant matrix and V be an "incomplete 
Vandermonde matrix," i.e., the 2Vj x TV matrix 

v=f; |. (2.6) 

Uf- 1 ... 

Then, the soliton potential is given by 

u(x) = -2d 2 Xi \ogr(x), (2.7) 
where the r-function can be expressed as 

t(x) =det(Ve^X»). (2.8) 

For the Jost and dual Jost solutions we have 

$(x,k) =e- lkxi - ] * 2x2 x(x,k), (2.9) 
V(x, k) = e lk *i+ k2 «»£( Xj k), (2.10) 



where 

X(*,k) = ^, (2.11) 

^, k ) = Z*^), (2 .i 2) 
t{x) 

with (Miwa shift) 

/ e K{x) \ 

(s,k) = dct(Ve K ^(K + ik)£>), r*(x, k) = det f V — — 2? j , (2.13) 



and 



« + ik = diag{«n + ik}^L x . (2.14) 

In order to study the properties of the potential and the Jost solutions, it is convenient 
to use an explicit representation for the determinants. By using the Binet-Cauchy formula 
for the determinant of a product of matrices we get 

, A/" N b 
T W = JH E n{ni})V{{ ni })Y[e K ^\ (2.15) 
b ' {m}=i i=i 

x(x ' k)= AM^ E mmi})V({ mi })l[(K mi +ik)e K ^, (2.16) 

«*. k) = ^ ^ £ i WW) f[ (2.17) 
where we used notation 

/I ... 1 \ 

^({n,}) = det : : =11 ( 2 - 18 ) 

V K^- 1 ... K^- 1 / 1<*<J<*» 
\ "1 ™JV b / 

: : (2.19) 

for the maximal minors of matrices V and 2? and where 

{mj} = {mi,...,mN h }, {rii} = {n%, . . . ,riN b } (2.20) 

stand for non ordered sets of Nb indices from the interval 1, . . . , AT. 

Notice that the only x-dependent terms in (12.15[) . (|2.16l) . and (|2.17l) are exponents of 
sums of linear functions (|2.4[) . Correspondingly, the asymptotic behavior of the function 
t(x) and of the potential has a sectorial structure on the x-plane. In order to specify these 
sectors at x — > oo we introduce the ray directions 

xi + (K n + K n+Nb )x 2 bounded n — \ J\f (221) 

{n n +N b - K n )x 2 -)■ -OO. ' )•••: J 

where we assume that the indices are defined mod TV, so that thanks to (|2.2p . say, n + Nf, — 
n — N a for n > N a . Thus there are N a rays in the direction x 2 — >• — oo and Nb rays in the 
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direction X2 — ► +00. The sector a n is swept out by rotating anticlockwise the ray r n up 
to the ray r„ + i. These sectors are nonintersecting and cover the whole x-plane with the 
exception of rays. In we proved that the the leading exponents of t(x) when x — > 00 
are the exponents exp(^™ = ^ V6 1 Ki (x)) , each being the leading one in the corresponding a n 
sector of the x-plane. More exactly, if the coefficients 

z n = V(n n , n n+Nb ^)V{n, ...,n + N b -l), (2.22) 

are different from zero for all n = 1, . . . , Af (again with indices defined mod TV) the function 
t(x) has the following asymptotic behavior along rays and inside sectors: 

/n+N b -l \ 

x^oo: t(x) = (z n + z n+1 e K ^^- K ^ +o(l))expi Jl K l( x )h ( 2 - 23 ) 

/n+N b -l \ 

x^oo: t(x) = (z n + o(l))exp( ^ Ki{x) . (2.24) 

Regularity of the potential u(x) on the x-plane is equivalent to the absence of zeroes of 
t(x). It is clear that it is enough to impose the condition that the matrix T> is Totally Non 
Negative (TNN), i.e., that 

T>(n 1 ,...,n Nb ) > 0, for all 1 < n x < . . . < n Nb < N. (2.25) 

However, sufficient conditions on the matrix T> for the regularity of the potential are un- 
known. On the other side, from (|2 . 23[) and (|2 . 24[) it follows directly that it is sufficient to 
require that 

z n > (2.26) 

for having nonsingular asymptotics of the potential. 

We also mention that the functions x(a;,k) and £(x,k) have bounded asymptotics on 
the x-planes because the x-dependent exponents enter in denominators and numerators of 
expressions (I2.16[) and (|2 . 1 7[) with coefficients proportional to V({rii]). This means that the 
leading asymptotic behavior of the denominators of the functions x( x i k) and £(2;, k) on the 
x-plane is not weaker then the behavior of their numerators. For more details see |12j . the 
review papers [El [14] , and [15j [ITJ [16] , where the same notations have been used. 

We need in the following also the values x( x ? i^n) of xi x i k) at k = k n and the residues 
i n {x) of £(x, k) at k = iK„. From (EU, (I2TTO1) and (j2~!B . (I2TT7)) we have 

x ( X) i Kn ) = L> Yl ^({™j)n{^},«)n^ (a;) , (2.27) 

b- T \ x ) { m .} =1 i = i 

Ux) = ~NTT-) £ mni})J^S njn (-iy-W( ni ,...,n],...,n Nh )l[e K ^\ (2.28) 

1 b - T[ - X > {„,}=! j=l 1 = 1 

where {{rrii},n} = {rn\, . . . , rriN b , n}, hat over rij denotes that this index is omitted and 
where the 8 njn Kronecker symbol in the r.h.s. of the last formula is due to the fact that the 
residues of the terms in the sum are nonzero only when some rij = n. 

Taking into account the analyticity properties of x( x > k) and £(x,k) in (|2.16[) . (I2.17[) 
their product can be written in terms of the values X (x,iK n ) and £, n ( x ) as follows 



x(x,iK n )^ n (x') 

n=l 



X (x,mx',k) = 1 + T ^/^"^ (2.29) 
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which also will be useful in the following. In [16. we demonstrated that the Jost solutions 
obey the Hirota bilinear identity 

M 

^ / $(x,iK n )* n {x')=0 ) (2.30) 

n=l 

where in analogy to (I2.28|) ^ n (x) denotes the residue of ^(x,k). 



3 Green's function 

We want, now, to prove that 

Q(x, x',k) = _ s § n ( x 2-4) f ds & ^ s 2 _ k | )(a , 2 _ s + i k a )*(x', a + i ka) + 



2tt 

+ i6(x' 2 - x 2 )^26(k<3 -K n )${x,iK n )$ n (x'), (3.1) 

n=l 

where 9 denotes the step function of its argument, satisfies (II ,3|) and that G(x, x' , k) defined 
in (|1.4I) is bounded with respect to the variables x,x' el 2 and k € C and has finite limits 
at space infinity. 

Notice that convergence of the integral in (|3.1j) follows directly from boundedness of 
functions x(x, k) and £(x' , k). Applying the heat operator (|1.2[) to 0(x, x' , k) in (|3.1|) we get 

C(x,d x )G{x,x' ,k) = ^lZfM f ds $fa g + jk Q ,)*(a; / ! s + ik9)+ 

2n J 

AT 

+ iS{x 2 - 4) 0( k9 i«n)*n(a:'). (3-2) 

n=l 

The integral, after inserting (|2.29[) in it, can be explicitly computed getting 
* 1 " 2 ' fds §{x, s + iksj)*(a;' ', s + ik®) = 8{x- x')- 

- id(x 2 - x' 2 ) ^ #( k o -n n )®(x, in n )^ n {x') 

71=1 

+ i<5(x 2 - x' 2 )6(x\ - H) J] $0, *Kn)*n(»'). (3-3) 



2tt 



Then, thanks to (|2.30p . equation f|l .3[) is proved. 

In order to prove boundedness of G(x, x' , k) we write 

G(ar, at', k) = G c (z, a;', k) + G d (z, i', k), (3.4) 

where the "continuous" part G c (x, x' ,k) of the Green's function equals 

G c (x,x',k) = - Sgn{x *~ X ' 2) Jdse^-k^ix^x^e^-^-^+^^-^x 

x e (k « -^'-"^(as, s + i k^(x', s + i k 3 ). (3.5) 
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and its "discrete" part Gd(x, x' , k) is given by 

G d (x, x', k) = i9(x' 2 - X2 yK^-<)+^ ^2 0(k 9 -Kn)$(x, iK n )* n (x'). (3.6) 

n=l 

While the boundedness of G c (x, x',\t) in (I3.5[) for cc, a;', and k (including limits at infin- 
ity) follows directly from boundedness the product %(a;, k)£(a/, k) at infinity, properties of 
Gd{x, x', k) in (|3.6I) deserves a more accurate investigation. 

Let us consider the sum in the r.h.s. of ([376]) . Inserting in it (1237)) and (|2T28|) . thanks to 
antisymmetry of minors of matrices V and V (see (f2~T8]) and (|2?l9|) ). after summing over n, 
we get 

^ 0(k 3 -K n )$(x, iK n )* n (x') = 
ri=l 



7T E E ©({m i })2?({n i })d(ka-/ Sri 



Nb\(Nb - 1)\t(x)t(x') w=1{ni}=1 

X ^({^},«ArJ^(«l,---,«Ar b -l)>< 

-1 \ 

(3.7) 



/JVi, ^b-l \ 

exp K ™i s ) + K «« b ( x ) + E 

\i=i i=i / 



We recall, now, that the maximal minors of a matrix satisfy the Phicker relation, i.e., 
for any subsets {mj and {n^} of indices running from 1 to TV and arbitrary j 6 {1, . . . , Nb} 

V{{mi})V{{ni}) = 

N„ 

= y^P(mi, . . . ,m s -i,nj, m s+ i . . . ,mN b )T>(nx, . . . ,nj-i,m a ,nj+i,. . . ,n Nb ). (3.8) 

s=l 

Inserting this equality with j = Nb in the r.h.s. of (I3.7[) we exchange m s njq b for s — 
1, . . . , Nb- Notice that under this transformation the first Vandermonde determinant changes 
sign, while the second Vandermonde determinant is unchanged, as well as the exponent. 
Thus 

A/" 

#(ka -K n )$(z, in n )^ n {x') = 

n=l 

N b Jif A/" 

wT^yE E E V{{ mi })V{{n % })* 



N b \{N b - 1)\t(x)t(x>) ^ ^ ^ ■ 

s=l {rrii} = l {n;}=l 

x 6>(kcj -K ma )V {{mi} ,n Nb )V '(ni, . . . , rc^-Ox 

/ N b N b -1 \ 

x exp ^ (x) + K nNb (x) + K n> ■ 
\i=i i=i ) 



(3.9) 



Exchanging now m s O m N we get, summing over s, iV;, equal terms. Finally, we multiply 
(|3.7p by TV;,, sum up with (13.91) . divide this sum by N b + 1 and insert the result into (13.61) . 
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This gives us 
G d (x,x',k) = 



^ '-• ) - J2 E »(W)»({«i})^ mw )x 



- Dim + iM.MxO {m ^ =1{n ^ =1 

x [6>(k 9 -« ) - 0(k 9 -K nN )]V({m,i},n N )V(n u ...,n N x )x 



AT, JV.,-1 



xcxp ^ # mi (x) + (x) + ^ (a;') + i k(a; 1 - x[) + k 2 (x 2 - x' 2 ) , (3.10) 
\i=i 1=1 / 

where 

k mn = k? R -(kcj -K m )(kcj -K n ), (3.11) 

and, for convenience, we introduced the factor 9{k m ^ n ), which identically equals 1 for kg in 
between K m and K n , since fc TO; „ is positive inside this interval. Let us also denote 

Zmn = Xl + (n m + K n )x 2 , (3.12) 

and analogously for z' and x', and let us notice that 

fl(kcj - « m ) - 0(kcj -« n ) = sgn(z m „ - z' mn )x 

x [0((k 9 -K m )(z mn - O) - (( k a - K «)(^ m „ - O)]' ( 3 - 13 ) 

and that 

— x^) + k (x 2 — x 2 ) + K m (x) = k mn (x 2 — x 2 ) — (kg —n m )(z mn — z mn )-{- 

+ ik SR [x 1 - x\ + 2k a (x 2 - x' 2 )\ +K m (x'). (3.14) 

Then, we can decompose Gd(x, x', k) in the sum 

G d (x, x', k) = Gi(x, x', k) + G 2 (x, x', k), (3.15) 

with 

Gi(x,x',k) = 



giks R [x 1 -a;' 1 +2k3(2; 2 — x' 2 )\ ^ 

((JVft - l)!) 2 (iVb + l) E E S S n (~m.v -»'v„»v,, 
{mj}=l {n^}=l 



X @{km Nb n Nb (%2 x 2))^((k3 K n Nb )(Zm Nb n Nb z rn Nb n Nb )) x 

x e fcm^»^(^-^)-(ka -^)(*» W6 -*m Wfc „ Wfc V({m i },n JV J^(n 1 , . . . ,njv fc _i)x 
D(M) exp ^ X m! (x)^J V({m}) exp ^ /f n , (x')^J 



(x) r(x') ' (3A6) 
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and 



G 2 (a;,x',k) 



iU 3t [x 1 -x' 1 +2\t a (x 2 -x' 2 )] 



Ml JV JV 



((A 6 - 1)!) 2 (A 6 + 1) ^ ^ s g n U»« 6 «« 6 'm^n^, 
{m l }-l {ni} — 1 

X 9(km Nb n Nb (x 2 — X2))d[(k.Q ~Hm Nb ){Zm Nb n Nb ~ z m Nb n Nb )) X 
x e k ™N b KN b ( x 2- x 2)-(ka - K m Nb )(Zm Nb n Nb -z' mN ^ nN ^ ^ 

/N b -1 \ 

2?({rrii})exp( K mi (a;) + K nNb (x) J 
x V{{mi},n Nb ) 



f = l 



r(x) 

exp ( jr K ni (a/) + (,')) 
xV(n 1 ,...,n JVi _i) ^ — ; (3.17) 

T(X ) 

The exponents in the third lines of the r.h.s. of (|3.16[) and ()3.17p are decaying or at least 
not growing thanks to the ^-functions in the second lines. So we have to check behavior 
with respect to x and x' of the last lines of these equations. About A3. 16)) the situation is 
trivial. The two exponents in the forth line has the same coefficient (minor of T>) as in t{x) 
and t(x') in the denominator and, therefore, they are bounded when x and x' are growing. 
Situation with (|3.17[) is more involved. Let us consider the term in the forth line. If the 
minor T>({rrii\) in the numerator is different from zero and "D(n%, ■ ■ ■ , njv 6 _i) ^ 0, then the 
same exponent as in the numerator is present in r(x) in the denominator and the ratio is 
bounded. But if T>(nx, . . . , riN b -i) = such exponent is not involved in t(x) and in the 
direction where it is the leading one (if such direction exists) the ratio is growing at large 
space. The same is valid for the term in the fifth line of (|3.17|) . Thus, in the case of a Totally 
Positive matrix V the boundedness of the Green's function G(x,a;',k) (|1.4[) is proved. If 
instead of TP we impose on the matrix V conditions (|2.26p , then all leading exponents are 
involved in t(x), as we mentioned in discussion of (12.23[) and (|2 . 24|) . Thus again, both ratios 
in (I3.17P have finite asymptotics at space infinity, but can be singular in the finite domain. 
To avoid this singularities it is enough to impose additionally that the matrix T> is TNN. 
But the case of a generic TNN matrix T> requires an additional study. 

Boundedness of the Green's function G(x, x 1 , k) with respect to the variable k also follows 
from (|3.5p . p. 161) and (|3.17p . In all points k = i/ci, . . . this function is discontinuous. 
As it was shown in j9_, in the case N a = Nb — 1, in order to give the detailed behavior of the 
Green's function at the points of discontinuity and in order to determine possible additional 
Green's functions useful in defining spectral data and studying their properties one needs 
a more generic object: the extended resolvent of the heat operator. Here, in the following 
section, we introduce the extended resolvent, leaving for a forthcoming article the study 
of its properties (much more involved than in the case of the nonstationary Schrodingher 
operator) and its use in building the inverse scattering theory. 

4 Extended resolvent 

Let us recall that the resolvent M(q) of the extended heat operator L(q) which has kernel 
L(x,x';q) =£{x,d x + q)S{x-x'), q - {q u q 2 ) £ R 2 , (4.1) 
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is defined as the operator with kernel M(x, x'; q), which is the left and right inverse of L(q). 

In [16] we showed that for a pure soliton potential the operator L(q) has left annihilators 
if N a < Nf, for q belonging to the union of some polygons included in the polygon inscribed 
in the parabola q2 — qf of the g-plane, with vertices at the points (k„, re„), n — 1, ... ,7V 
and, analogously, we showed that, for N a > Nf,, the operator L(q) has right annihilators 
in an analogous region of the q-plane. Therefore, we need to find an operator M(q) with 
bounded kernel M(x,x';q) defined for any q £ M 2 , which, for N a < Nf,, is right inverse of 
L(q) for any q, but left inverse only for q not belonging to the union of polygons quoted 
above, and which, for N a > Nb, is left inverse of L(q) for any q, but right inverse only for q 
not belonging to an analogous union of polygons. Here we give only the result and we leave 
the lengthy and involved deduction to a following paper. 

The searched operator M(q) has kernel 

M(x,x';q) = M c (x,x';q) + M d (x,x';q), (4.2) 

where 

e -q(x-x') l- 

M c (x,x';q) = -sgn(x 2 - x' 2 ) — d Pl 8((q 2 + p\ - qj)(x 2 ~x' 2 ))x 

x<f>(x;pi+iqi)V(x';pi+iqi), (4.3) 



and 



{ rn } — 1 {n^} — 1 

x [6(q! -K mNb ) -0(q% - K nNb )}9(q mNb n Nb (x 2 ~x' 2 ))x 
x V({m,i},n Nb )V(ni, . . . ,n Nb -i)x 

/ N b N b -1 \ 

x exp I ^ ^ TO1 (x) + K nNb (x) + J2 K m i x ') . ( 4 - 4 ) 
\i=i i=i / 



with 



fan = q.i - (K m + K n )qi + K m K n . (4.5) 
One can check directly that the Green's function Q(x, x', k) can be obtained via the reduction 

G(x,x',\c) = e i ^ x i-< +2 ^ x *- x '^M(x,x';kz,kl-kl). (4.6) 

Notice that the points (qi = kcj,^ = k Q — k^) lies outside the parabola q 2 = qf in accor- 
dance with the fact that in this region the operator L(q) has right and left inverse. On the 
parabola the only points, where L(q), for N a ^ Nb, has no inverse are the points k = k,„ 
n = 1, ... ,7V and, therefore, as expected, the Green's function G(x, x', k), being a reduction 
of the resolvent, is discontinuous just at these points. 
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